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I. INTRODUCTION 

During recent years, a lot of work has been devoted to 
test the Standard Model (SM) of elementary particles. 
The best tested are interactions between gauge bosons 
and matter, and in this sector the SM description turns 
out to be very accurate. Other sectors, however, has 
been checked to much less degree. Among them are self- 
interactions of gauge bosons as well as the Higgs sector, 
which plays important role for completeness of the model 
and in many aspects of symmetry breaking. Also many 
shortcomings of SM, like the big number of free parame- 
ters, unresolved question of mass hierarchy, the problem 
of massive neutrinos and their oscillations, may call for 
more desirable description of Nature. 

In the matter of fact, a number of various models 
reaching beyond SM's orthodoxy were proposed. One of 
the most promising candidate is the supersymmetric ex- 
tention of SM called Minimal Supersymmetric Standard 
Model (MSSM). It is based on the concept of supersym- 
metry (SUSY) and, despite the lack of direct experimen- 
tal evidence at the moment, is supported by many theo- 
retical arguments accompanied with hope, that SUSY is 
the relevant description of our world above 1 TeV scale. 
One of the main facts supporting MSSM is that incorpo- 
rating SUSY in SM causes all the gauge couplings unify 
at some scale ttigut ~ 10 16 GeV. As is well known, ex- 
trapolations of data from the LEP measurements suggest 
such behaviour. However, SUSY particles have not been 
observed in experiments, so supersymmetry has to be 
broken in low energy regime. The issue how this break- 
ing is realised, is the least understood question of the 
theory. 

The most widely studied version of SUSY conserves 
the so-called ii-parity. The ii-parity is a multiplicativ 
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quantity defined as R = (_i) 2 -S+3B+£ ; where B and L 
are the baryon and lepton numbers, and S is the spin 
of corresponding particle. As a consequence, processes 
which do violate lepton or baryon number are strictly for- 
bidden unless the symmetry is broken. Moreover, SUSY 
particles are pair produced and the lightest one is stable. 

The origin of ii-parity conservation is not based on 
any fundamental principle, so this property of MSSM is 
an ad hoc hypothesis and therefore some extentions of 
the model, allowing the violation of ii-parity, were dis- 
cussed in literature. These modifications can be classified 
either as explicitely ii-parity broken MSSM (j?MSSM) 
approaches [ij or as formalisms with spontaneous break- 
ing of this symmetry 2J. In the first class of models the 
R-violating interactions are consistent with both gauge 
invariance and SUSY Q while the second ones provide 
the simplest way for ii-parity violating effects conserv- 
ing at low energy the baryon number [4(. The explicit 
ii-parity breaking leads to well defined phenomenologi- 
cal consequences, but due to large number of free param- 
eters involved, such theory has only marginal predictive 
power. In contrast, the spontaneous breaking has many 
virtues added, like the important possibility of dynamical 
origin of the ii-parity breaking [4( . 

Theories of gauge mediated supersymmetry breaking 
(GMSB) belong to the second kind of approaches and 
have recently attracted a great deal of attention. They 
are highly predictive, offer a natural solution to the 
flavour problem and contain much less free parameters 
compared to MSSM with SUSY breaking mediated by 
gravitational interaction 0, 0, 0, S 0, E3, E3, E2| In 
GMSB models supersymmetry breaking is transmitted 
to the superpartners of quarks, leptons and gauge bosons 
via the usual ST/ (3) x SU(2) x U(l) gauge interactions 
and occurs at the scale Msusy ~ 10 5 GeV. Gauginos 
and sfermions acquire their masses through interactions 
with the messenger sector at one- and two-loop levels 
respectively, resulting in different phenomenology of the 
low-energy world from the MSSM one. In these mod- 
els flavour-diagonal sfermions mass matrices are induced 
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in a rather low energy scale, and therefore they supply 
us with a very natural mechanism of suppressing unob- 
served in experiments flavour changing neutral currents 
(FCNC). Moreover, since the soft masses arise as gauge 
charges squared, the sizeable hierarchy proportional to 
the gauge quantum numbers appears among the super- 
partner masses. In this light, recently renewed interest 
in GMSB |<llic| is understandable. 

The i?-parity in MSSM can be explicitely violated by 
the presence of bilinear ^| and trilinear [T^j terms in the 
superpotential. The trilinear terms lead to lepton num- 
ber and flavour violation, while the bilinear terms gen- 
erate non-zero vacuum expectation values for sneutrino 
fields (£;), causing neutrino-neutralino and electron- 
chargino mixing. Thus, approaches with lepton number 
violation can describe some low-energetic exotic nuclear 
processes like the neutrinoless double beta decay (0vf](3) , 
known to be very sensitive to some of the i?-parity vio- 
lating interactions |l5j]. Using experimental data about 
these processes, e.g. bounds on the half-life of neutrino- 
less double beta decay, one can establish str inge nt limits 
on the ^-parity breaking SUSY [Tl 111 ITl ITa liH] . 

Supersymmetric models with i?-parity non- 
conservation have been extensively discussed in the 
last decade (see e.g . fl3l U^\ - and were also used for 
the study of [iH~EH~EE 13 13 13 . The older 

calculations were concentrated on the conventional 
two-nucleon mode of 0v(3(3, in which direct interaction 
between quarks of the two decaying neutrons causes 
the process 0, 0. Recently the dominance of pion 
exchange mode based on the double-pion exchange 
between the decaying neutrons over the two-nucleon one 
was proved [19t |2CI l2l| . 

Motivated by the forementioned features of GMSB 
models, in this paper we study the i?-parity breaking 
phenomenology of MSSM, and use the neutrinoless dou- 
ble beta decay for deducing limits on some non-standard 
physics parameters. In the previous studies such esti- 
mates were performed in the framework of |?MSSM with 
supergravity mediated SUSY breaking by means of GUT 
constraints 0, 0] or additional assumptions relating 
sfermions and gauginos masses 0, [l^. We will show 
that one can find quantitatively new constraints [l8j 
within GMSB models. In this paper we study this prob- 
lem using up-to-date experimental data from CLEO col- 
laboration [22] for 0z//3/3. As previously, we limit our at- 
tention to the trilinear terms only, leaving complete treat- 
ment of bilinear and trilinear i?-parity violating terms in 
GMSB for subsequent paper. 

For reliable extraction of the limits on i?-parity break- 
ing constant A' m from the best presently available exper- 
imental lower limit on the half-life of 0vf3{3, it is neces- 
sary to determine other SUSY parameters, e.g. masses 
of SUSY particles, within a proper SUSY scenario, and 
to evaluate corresponding nuclear matrix elements. Be- 
cause at present the Renormalize d Quasip article Random 
Phase Approximation (RQRPA) [23, [24j , which takes to 
some extent the Pauli exclusion principle into account, is 



the main method commonly used in calculations of the 
0vf3(3 nuclear matrix elements [2l| , we used it also in our 
work. 

Our article is organized as follows. In Section II 
the necessary theory is developed. We also discuss to 
some extent the gauge mediated supersymmetry breaking 
mechanism of the neutrinoless double beta decay. Section 
III contains the results and analysis of constraints im- 
posed on supersymmetric parametres by non-observation 
of Qvj3(3 in Germanium 76 Ge, for which the best exper- 
imental limit on the half-life is known. In this part we 
also demonstrate differences between the neutralino and 
gluino mechanisms in the neutrinoless double beta decay. 
Finally, summary and concluding remarks can be found 
in Section IV. 



II. FORMALISM 

A. 7?-parity violation in MSSM 

In this section we briefly outline main features of 
MSSM and its violation machanism. Both in the super- 
gravity and in GMSB, the i?-parity can be explicitely 
violated by the bilinear [l3l | and trilinear |l4j terms 
incorporated into the superpotential. Bilinear terms 
generate non-zero vacuum expectation values for the 
sneutrino fields (vl), causing neutrino-neutralino and 
electron-chargino mixing. Trilinear terms lead to the 
lepton number and flavour violation. Above features 
make |?MSSM models appropriate for the description of 
OvP/3. Because this process is known to be very sensi- 
tive to supersymmetric and i?-parity breaking parame- 
ters, data from the nowadays double beta experiments 
allow to establish stringent limits on ^?MSSM physics 

[3l3l3l3l3l3l2illll 

The complete superpotential W of the model can be 
written in the form 

W = W + Wp, (I) 

where 

Wo = hf j Q i H u u c j + hP j Q i H d d c j 

+hf j L i H d e c j + fiHdHu (2) 

and 

Wft = A,,,. I., I., ill + X' ijk ,<!'., 

+\'l jk u c l d c J d r k + H L J H u (3) 

are the i?-parity conserving and i?-parity breaking parts, 
respectively. Here Q and L denote the quark and lep- 
ton SU(2) doublet superfields, ii c , d c , and e c the corre- 
sponding SU(2) singlets and H u , Hd are the Higgs su- 
perfields. In the i?-parity breaking part ©, the two first 
terms are lepton number violating while the third violates 
the baryon number conservation. The presence of these 
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terms simultaneously would cause unsuppressed proton 
decay and therefore we follow the usual way and simply 
set Xijk = \'(j k = in order to avoid such possibility. 

In the low energy world supersymmetry is obviously 
broken and usualy one supplies the theory with the "soft" 
breaking terms, being another source of i?-parity viola- 
tion: 

= (AfjQiHuUj + APjQiHddj + Af j L i H d e c J + h.c. 
+ Bfi (H d H u + h.c.) + m 2 Hd \H d \ 2 + m 2 Hn \H U \ 2 
+ ml\L\ 2 + m| |e c | 2 + m 2 ~\Q\ 2 + mtju c \ 2 + m 2 ~ \d c \ 2 

+ Qm^bVs + ^M 2 $frl>w + \m~ a i> a g ^ a g + h.c. 



(4) 



and 



■■soft 



\ijkLiLjU c k + X' ijk LiQjd k 



ijk 

\jk a i u j u k T P^jLjH u 



Mij Lj Hd ■ (5) 

Here, fields with tilde denote the scalar partners of quark 
and lepton fields, while ipi are the spin-^ partners of 
gauge bosons. 

To describe Qv[3f3 process within supersymmetric mod- 
els one needs an explicit form of the appropriate La- 
grangian. It can be obtained using the standard proce- 
dure of extracting Lagrangian from superpotential Wn. 
After some computation one gets 



(uL,d L ) 
+ (ei, v L )d R 



-d% 



(6) 



+ {u Ll d L )d t 



h.c. 



Applying the formalism described in details in e.g. [iH 
Il9j . one ends up with the effective Lagrangian: 



VPS JpsJpS — -£)T Jt Jthv 



(7) 



where the color singlet hadronic currents are Jps — 
u a (l + j5)d a , Jy" = ^"(^(l + 75)c? a , with a being the 
color index and — (z/2)[7 M , 7"]. The effective lep- 
ton number violating parameters rjps and tjt in Eq. (|7|) 
accumulate fundamental parameters of MSSM and their 
explicit forms, obtained with a proper treatment of the 
colour currents in the Lagrangian, can be found in Ref. 
|19|. These parameters are rather complicated functions 
of supersymmetric masses and coupling constants which 
are, in general, free quantities limited by experimental 
data or theoretical considerations only. In the next sec- 
tion we describe our procedure how to obtain values of 
most of them in the GMSB MSSM model. 



B. GMSB MSSM and procedure for finding 
supersymmetric parameters 

Supersymmetry breaking in GMSB models occurs in 
the so-called hidden (or secluded) sector. It is a well 
known fact, that the detailed structure of this sector does 
not change the phenomenology of low energy world. In 
our approach we assumed that the secluded sector con- 
sists of a gauge singlet superfield S, whose lowest S and 
F components acquire vacuum expectation values (vev). 

Supersymmetry breaking is communicated to the visi- 
ble world via the so-called messenger sector. The inter- 
action among superfields of the secluded and messenger 
sectors is described by the superpotential 



W = X t S^ 



(8) 



where $i and <&i denote appropriate messenger super- 
fields. Because of non-zero vev of lowest S and F compo- 
nents of superfield S, fermionic components of messenger 
superfields gain Dirac masses M, = XiS and determine 
in this way the messenger scale M. Simultaneously mass 
matrices of their scalar superpartners 



\XiS\ 2 
X*F* 



XiF 

\x z s\ 2 



(9) 



have eigenvalues \XiS\ 2 ± \\iF\. 

It is easy to see that vev of S generates masses for 
fermionic and bosonic components of messenger super- 
fields, while vev of F destroys degeneration of these 
masses, which results in supersymmetry breaking. Defin- 
ing Fi = XiF one can introduce a new parameter Aj = 
Fi/S measuring the fermion-boson mass splitting: 



nif = Mi . 



A, 



(10) 



Parameter A and messenger scale M are in the following 
treated as free parameters of the model. 

Messenger superfields transmit SUSY breaking to the 
visible sector. It is realized through loops containing in- 
sertions of S and results in gaugino and scalar masses at 
M scale: 



M- X .(M) = fc^A G , 



4k 



(11) 



7 



\{M) = 2Y J C!k i (^]\% (12, 



4tt 



where i = 1,2,3 is the gauge group index, and 



N 



(13) 
(14) 
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with k being the flavour index. In Eqs. (|13|l and i|14|) 
is the doubled Dynkin index of the messenger su- 

perfield representation with flavour k. Coefficients C{ 
are the quadratic Casimir operators of sfcrmions. For 
d-dimensional representation of SU(d) their eigenval- 
ues are C = (d 2 — l)/2d. In the case of U(l) group 
C = Y 2 = (Q — T 3 ) 2 . It follows that coefficients ki are 
equal to 5/3, 1 and 1, for SU(3), SU(2), and U(l) re- 
spectively. The normalization here is conventional and 
assures that all kioti meet at the GUT scale. Finally, the 
functions / and g have the following forms: 

g(x) = -5 [(1 + x) log(l + a;)] + (x -> -x), (15) 



— Lin , 
2 VI 



log(l + a:) - 2Li 2 
2x 



+ (x 



l + x 
x). 



(16) 



The minimal model of GMSB considered in this paper 
contains only one messenger field flavour. Thus, dropping 
flavour indices, one can write 



M- Xi (M) = Nk 



mj(M) 



2N 



aj(M) 
4-rr 



, M 1 



(17) 



E?=i C(ki (^) 2 A 2 / (£) . 1, 



(18) 



where C[ = Y 2 , C{ = 3/4 for SU(2) L doublets and for 

singlets, C3 is equal to 4/3 for SU(3)c triplets and for 
singlets. In 118|l 1 denotes the unit matrix in generation 
space and guarantees the lack of flavour mixing in soft 
breaking mass matrices at messenger scale. N, the so- 
called generation index, is given by N = Ylj^i n-i, where 
N g means the total number of generations. In this paper 
we study two cases: 

1. a single flavour of 5 + 5 representation of SU(5), 
with SU{2) L doublets [l and I), and SU(3) triplets 
(q and q); 

2. a single flavour of both representations 5 + 5 and 
10 + 10 of SU(5) group. 

In case 1. AT is equal to 1, while in case 2. N = 1 + 3 = 4, 
because for 10 + 10 representation of SU(5) the doubled 
Dynkin index is 3. 



C. Renormalization Group Equations and 
parameters determination 

The evolution of all running parameters is realized us- 
ing Renormalization Group Equations (RGE). The for- 
mulae (|17f) and (|18fl may therefore serve as boundary 
conditions for evolution of soft parameters at the elec- 
troweak scale. Our procedure resulting in low energy 
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FIG. 1: Sample RGE evolution of sparticle masses in GMSB 
MSSM and their dependence on A. 



spectrum of SUGRA and GMSB MSSM models and its 
application to the description of 0v(3(3 decay can be found 
in our previous papers 0] , so here we only sketch its 
most important features. 

At the beginning, one evolves all gauge and Yukawa 
couplings for three generations up to the messenger scale 
M. We use the one-loop Standard Model RGE [2<J below 
the mass threshold, where SUSY particles start to con- 
tribute, and MSSM RGE ^ above that scale. We admit 
not to use the full set of RGE appropriate for the^?MSSM 
model [2^, [2?| . The influence of i?-parity breaking con- 
stants on other quantities running from the messenger 
to the electroweak scale is marginal due to the smallness 
of A's. In our case the two-loop corrections can also be 
safely neglected (for a discussion of this problem see |30()- 
Initially, scale Msusy is taken to be equal to 1 TeV, 
but it is dynamically modified during running of relevant 
masses. In the next step we construct the gaugino and 
sfermion soft mass matrices using Eqs. <|17[) and l|18|) . and 
perform RGE evolution of all the quantities back to mz 
scale. During this run, tt? h reaches negative value caus- 
ing dynamical electroweak symmetry breaking (EWSB). 
It is well known, that proper treatment of this mechanism 
needs minimizing of the full one-loop Higgs effective po- 
tential 31 . On the other hand, appropriate corrections 
contain functions of particle mass eigenstates generated 
by EWSB mixing. Thus, as the first approximation, we 
minimize the tree-level Higgs potential parameters \x and 
which are crucial for further analysis. 

Having all needed mass parameters at electroweak 
scale, one can evolve all other quantities to some scale 
Qmin, which is optimal for minimization of the one-loop 
corrected Higgs potential. At this scale, defined as the 
geometric mean of stop masses, minimization procedure 
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FIG. 2: Constraints on GMSB parameter space. 



results in /i and B[i values. Next, all the quantities are 
running back to mz scale. Iterating this procedure one 
obtains stable values of fj, and B[i and then low energy 
spectrum for the considered model. Only four parame- 
ters: A, M, tan/3 = v u /vd and sgn(/i) remain free. The 
quantities v u and Vd are vev's of H u and H^, respectively. 

In Fig. H a sample evolution of sparticle masses ver- 
sus the A parameter is shown. Other parameters were 
tan/3 = 3, M = 500 TeV, N = 1. One sees that the 
masses of squark and gluino depend heavily on A, while 
in the case of selectron and neutralino this dependence is 
much weaker. 



D. Restrictions on low energy spectrum 

To impose restrictions coming from the present theo- 
retical assumptions and phenomenological data on the re- 
sulting spectrum is a non-trivial problem. We would like 
to obtain limits on physics beyond SM induced by Oi//3/3 
experiments, consistent with constraints coming from: 

1. finite values of Yukawa couplings at the GUT scale; 

2. proper treatment of electroweak symmetry break- 
ing; 

3. requirement of physically acceptable mass eigenval- 
ues at low energies; 

4. flavour changing neutral currents (FCNC) phe- 
nomenology. 

Below we will briefly discuss these sources of additional 
constraints. 

The first requirement comes from the RGE evolution 
procedure. It is well known that running of the Yukawa 
couplings is sensitive to initial (i.e. at the electroweak 
scale) values determined by tan (3. For very small tan /3 
(< 1.8) the top Yukawa coupling may "explode" before 
reaching the GUT scale. It follows from the fact, that 
Yto P (rnz) ~ 1/ sin/3. Similarly, other couplings Yj, and Y T 
"blow up" before the GUT scale for tan /3 > 50 because 
they are proportional to 1/ cos (3 at electroweak scale. 
Such behaviour of the Yukawa couplings limits the range 
of tan/3 to the interval 2-50. 

Another theoretical constraint is imposed by the 
EWSB mechanism. In order to obtain a stable min- 
imum of the scalar potential, the following conditions 
must hold: 



^Bf > (\ri 2 + m 2 Hu )(\»\ 2 + m 2 Hd ) 



2Bfi < 2\fx\ 



(19) 



They are always checked in our procedure during RGE 
running, and points which do not fulfill these conditions 
are rejected (see Fig. |2 points marked "EWSB"). Next 
restriction comes from the requirement of positive eigen- 
values of mass matrices squared at the electroweak scale, 
and allows to find combinations of free parameters pro- 
viding the negative (forbidden) eigenvalues marked in 
Fig. Elas "v.e.v.". 

The most interesting set of constraints has its source 
in the FCNC phenomenology. Such processes, strongly 
experimentaly suppressed, limit upper values of different 
entries of the sfermion mass matrices at low energies (cf. 
IHSHESOSI)- Here we consider the B -> X sl 
decay onl y. The effective Hamiltonian for this process 
reads |32ll35l| 



H c ft — — 



4Gf k*k ^ 

K^Ktb 



V2 



4=1 



a(n)Pi(v), (20) 
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where K is the quark mixing matrix (CKM matrix) 
and Pj are the relevant operators taken from Ref. [35| . 
Among the Wilson coefficients Ci(/i) two: C7 and C&, 
are the most important for the analysis of impact of the 
SM and MSSM interactions. (The leading order and the 
next-to-leading order analysis of these interactions were 
discussed in [32, H(j-) In order to costrain the low en- 
ergy spectrum of supersymmetric models using FCNC 
processes, it is a common practice to define the parame- 
ter J?7, which measures the extra (MSSM) contributions 
to the B — > X,7 decay: 



R 7 = 1 



C 7 ' 



(O)extra 



(mw) 



C) [m w ) 



(21) 



where the index (0) stands for the leading order Wil- 
son coefficients and the superscript extra indicates SUSY 
(charged Higgs, chargino, neutralino and gluino) contri- 



(Q)extra 



and C. 



(0)SM 



butions. Explicit expressions for C 7 
can be found e.g. in Ref. . Constraints on allowed 
values of R 7 are induced from the present experimental 
limits on the branching ratio BR(_B — ► X s j) measured 
by CLEO collaboration 22]: 

BR{B -» X sl ) = (3.21±0.43 stQt ±0.27 sast ) x 10~ 4 . (22) 

The theoretical dependence of BR(i3 — > X s j) on R? con- 
fronted with such experimental data allows to make the 
following estimate: 

- 6.6 < R 7 < -4.4 or 0.0 < R 7 < 1.3. (23) 

Using the above restriction, one can exclude certain val- 
ues of supersymmetric parameters, which result in the 
i?7 coefficient outside the allowed region In Fig. |3 

such points are marked as "b — > s + 7" . 

Looking on Fig. [21 one sees that the constraints de- 
duced from FCNC phenomenology are very sensitive to 
the sign of fj,. The same behaviour was also observed in 
SUGRA MSSM model (see e.g. Ref. 19]) and is mainly 
due to sensitivity of charged Higgs and chargino contri- 
butions to R 7 on sign of the \i parameter. 

The additional dependence of R 7 on both tan j3 and A 
parameters is shown for positive and negative fj, in Fig. 
|3J In the case fx > the parameter R 7 grows up for 
smaller values of tan [3 and behaves in opposite manner 
for fi < 0. Moreover, in the latter case R 7 is, in general, 
bigger which results in more stringent restrictions. More 
detailed analysis is presented in Fig. where the most 
important impacts to R 7 for different choices of tan /3, 
A, and sgn(/x) are explicitely shown. One can see that 
tan (3 and sgn(/z) do not influence the charged Higgs con- 
tribution significantly. Thus, a crucial point in the analy- 
sis becomes chargino contribution. Contrary to SUGRA 
MSSM case 0] the magnitude of chargino influence on 
i?7 is almost equal to the influence coming from charged 
Higgses. For positive values of the fi coupling constant, 
the chargino impact grows with increasing of tan /3, while 
for fi < case one observes opposite behaviour. In this 
light, behaviour of the surfaces shown in Fig. [21 becomes 
clear. 
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FIG. 3: R T parameter in GMSB MSSM for both signs of fi. 
The scan is performed over A and tan /3, with M — 500 TeV 
and N = 1. 



III. NEUTRINOLESS DOUBLE BETA DECAY 
AND LIMITS ON NON-STANDARD PHYSICS 



Restrictions imposed on the model by low energy spec- 
trum allow for reliable analysis of exotic nuclear pro- 
cesses, like the neutrinoless double beta decay, and then 
for deduction of additional constraints imposed on non- 
standard physics. In this paper we use experimental in- 
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FIG. 4: Contributions to c^' MbbM in GMSB MSSM comine from charged Hissscs and chareinos for both sisns of u. All 



8 



0.45 



i 1 r 

2n-mode 

Ji-mode 




100 150 200 250 300 350 400 
A [TeV] 



1.8 



T 



T 



T 



— i 1 r 

g-mechanism 

%-mechanism 




150 200 250 300 350 400 
A [TeV] 

FIG. 5: Contributions to A' ln coming from two-nucleon and 
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neutralino mechanism of SUSY breaking. All parameters like 
in Fig. ^except tan (3 = 3 and fi > 0. 



formation about non-observability of the Ov/3/3 decay in 
different nuclei to extract stringent limits on i?-parity 
breaking. 

The half-life of the process, taking into account all 
three possible t ype s of hadronization (2-nucleon, 1-pion 
and 2-pion plillflj) reads: 



/2J 



3 



i) T M.\ 
5 



2N 



Vps-Vt)M} n 



Vps M 



tN 



(24) 



In this equation M"~ N ,M 2 j: N and M* N are matrix el- 
ements for the 2N, lir and 27r channels, respectively. 



These matrix elements depend on non-standard physics 
parameters, involved in description of the neutrinoless 
double beta decay, and on nuclear structure details of 
decaying nuclei. (The explicit forms of elements (|24|l can 
be found, e.g., in [Hi!!.) Our procedure limits the 
number of free parameters to A, M, tan/3, sign(/i), and 
N only. As the loop diagrams with messenger fields do 
not affect the A-terms considerably, we can equal the 
common soft SUSY breaking parameter Aq to at the 
M-scale. 

Following well established procedure, the nuclear 
matrix elements in question were calculated within 
the proton neutron Rcnormalizcd Quasiparticle Random 
Phase Approximation (pn-RQRPA). This approach in- 
corporates the Pauli exclusion principle for fermion pairs 
|23ll24T| and is suitable for studies of nuclear structure as- 
pects of various double beta decay channels in open shell 
systems. Details of the method and its application to the 
double beta dec ay w ere presented in a number of articles 
(see, e.g., 00 El). 

Having both supersymmetric spectrum and nuclear 
matrix elements, one can extract from Eq. I|24|l con- 
straints on i?-parity breaking in GMSB MSSM using 
experimental information from non-observability of the 
neutrinoless double beta decay. Such approach is based 
on comparison of the theoretically obtained half-life for 
this process, as a function of some free non-standard pa- 
rameters, with the experimental upper limit for T1/2 in 
the given nucleus. 

We start with a presentation of constraints on A'n 
coming from different channels of Qv(3(3. In Ref. 
the problem of the pion mode has been discussed in de- 
tails. In Fig. [S] the importance of pion-exchange mode 
is clearly visible. The curve corresponding to the pion 
mode lies definitely below the line corresponding to the 
nucleon channel, so the pion mode imposes more strin- 
gent restrictions on the coupling constant. Also the role 
of various mechanisms leading to SUSY breaking are pre- 
sented. These data were calculated for 76 Ge nucleus, 
for which the best experimental limits, coming from the 
IGEX collaboration, are known (3^. One sees that the 
SUSY breaking mediated by neutralinos sets more severe 
limits on A' nl than the gluino mechanism. 

Further analysis is presented in Fig. |SJ We 
have included most of the nowadays known ex- 
perimental data (see and references therein). 
The interesting thing is, that the combina- 
tions A' m /[(TO 9 ~/100GeV) 2 (TOg/100GeV)5] and 
A' m /[(m g /100GeV) 2 (m x i/100GeV)i] remain nearly 
unchanged within wide range of A's. This allows us to 
estimate the lepton number violating constant: 
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FIG. 7: Limits on A' in for different structures of messeneger 
sector (see text for details) for planned 0vf3/3 experiments. All 
parameters as in Fig. Q3 



scenarios in the case of expected sensitivity of planned 
neutrinoless double beta decays. Two different messen- 
ger sector structures are here taken into account: the 
5 + 5 representation (N — 1) and both, 5 + 5 and 10 + 10, 
representations (N = 4) of 5/7(5}. We include parame- 
ters for three new experiments The GENIUS- 
MAJORANA-GEM project is expected to reach sensitiv- 
ity of T 1/2 - 2.3 • 10 28 y for 76 Ge. The MOON experiment 
has T 1/2 ~ 1.3 • 10 28 y and investigates the 100 Mo nu- 
clei, and the EXO-XMASS experiment will be sensitive 
to values of the half-life up to around T X f 2 ~ 2.2 • 10 28 y 
for 136 Xe. The relevant results are presented in Fig. [7| 
It is worth noting that for N = 4 the allowed values for 
the lepton number violating constant are much higher. 
The most promising results can be expected from the 
MOON project, which may set the best constraints on 
the i?-parity violating coupling constant. 



FIG. 6: Limits in GMSB MSSM on various combinations of 
X'm and masses of SUSY particles coming from experimental 
lower bounds on the half-life of 0v/3f3 decay in different nuclei. 
The corresponding nuclear matrix elements have been calcu- 
lated using pn-RQRPA method and the bag model. Other 
parameters as in Fig. |3] 



and 




<2.73-10~ 3 . (26) 



These results lower the allowed values in the first case by 
around 15% (cf. [3). 

We study also constraints coming from different GMSB 



IV. CONCLUSIONS 

We have presented an analysis of the current experi- 
mental state of neutrinoless double beta decay in the lan- 
guage of gauge mediated Minimal Supersymmetric Stan- 
dard Model. Combining theoretical, phenomenological, 
and experimental data we obtained a set of constraints on 
various non-standard parameters. In particular, new up- 
per limits for A' m /[(m,j/100GeV) 2 (mg/100GeV)2] and 

Ai n /[(mg/100GeV) 2 (m x i/100GeV)5] were extracted. 
A detailed discussion of the Wilson coefficients, the SUSY 
contributions to it and its dependence on the whole al- 
lowed range of tan/3 and A up to 400 TeV was pre- 
sented. Also some preliminary studies related to three 
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new planned §v(3(3 experiments were performed. 
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